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Abstract. A rainbow subgraph of an edge-coloured graph has all edges 
of distinct colours. A random d-regular graph with d even, and having 
edges coloured randomly with d/2 of each of n colours, has a rainbow 
Hamilton cycle with probability tending to 1 as n — > oo, provided d > 8. 



1. Introduction 

An edge-coloured graph is a rainbow if no colour appears more than once. 
We will study rainbow Hamilton cycles in edge-coloured graphs with n ver- 
tices where the number of colours available is also n; thus a rainbow Hamil- 
ton cycle uses each of the colours exactly once. 

We consider in this paper only (random) regular graphs and (random) 
colourings where each colour occurs the same number of times. If each colour 
occurs on q edges, we thus have qn edges, and hence the vertex degrees are 
2q. We use the standard notation G(n, d) for a uniformly chosen random 
d-regular graph on n given (labelled) vertices. We will only consider the case 
d = 2q even. (Hence there is no parity restriction on n.) Having sampled a 
random graph G(n, 2q), we then randomly colour its qn edges by n colours 
(1, . . . , re, say) with q edges of each colour, again choosing uniformly among 
all possibilities. We denote the resulting randomly coloured random graph 
byG c (n,2q). 

Our main result is the following on randomly coloured random regular 
graphs. (For some related results on the random graph G(n, m), see Cooper 
and Frieze [3].) We say that an event holds with high probability (whp), if 
it holds with probability tending to 1 as re — ► oo. (All unspecified limits in 
this paper are for n — > oo.) 

Theorem 1.1. Consider the randomly coloured random 2q-regular graph 
G c (n,2q), with n colours and q edges of each colour. Then, whp, there 
exists a rainbow Hamilton cycle if q>4, and not if q < 3. 
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Recall that it was shown by Robinson and Wormald [16, 17] that G(n, d) 
whp contains a Hamilton cycle as soon as d > 3. In our setting, when 
d = 2q has to be even, we thus whp have Hamilton cycles, ignoring the 
colouring, when d > 4, but rainbow Hamilton cycles only when d > 8. It 
is nevertheless remarkable that whp some Hamilton cycle manages to pick 
up an edge of each colour in a random 8-regular graph, when there are only 
four edges of each colour to choose from. 

Remark 1.2. In a similar direction, Robinson and Wormald [18] showed 
that a random 3-regular graph with o{y/n) randomly specified edges whp 
has a Hamilton cycle passing through all the specified edges (and, more- 
over, in randomly prespecified directions). It has further been shown by 
Kim and Wormald [14] that a random 2q-regular graph whp has an edge- 
decomposition into q Hamilton cycles, provided q > 2. 

It is natural to ask whether, similarly, a randomly coloured 2g-regular 
graph with n colours and q edges of each colour, as above, whp has an edge- 
decomposition into q rainbow Hamilton cycles. By computing the expected 
number of such decompositions (similarly to the proof of Lemma 3.2 below), 
it is easily seen that this is whp false when q < 4. We leave the case q > 5 
(when the expected number tends to infinity) as an open problem. 

The proof of Theorem 1.1 is based on the small subgraph conditioning 
method introduced by Robinson and Wormald [16, 17], and further devel- 
oped in [12], [15], [19] and [13, Chapter 9]. However, for this problem we 
have to consider the colourings of the small subgraphs too, see Section 3. 

Acknowledgements. This problem was suggested by Alan Frieze during 
the Conference on Random Structures and Algorithms at Emory University, 
Atlanta, 1995; he was originally intended as a coauthor but later declined 
this. Consequently we are pleased to be able finally to dedicate this work to 
him, marking 10 years since its beginnings, in which he was involved, and 
60 since his. 

We also acknowledge the assistance of the Maple algebraic manipulation 
package for the variance calculations in Section 6. Although the proof we 
found can be verified by hand, Maple was instrumental in finding that proof. 

2. MULTIGRAPHS, A BIPARTITE GRAPH, AND TRAFFIC RULES 

As usual in the study of random regular graphs (with small degree), it 
is convenient to extend the study to multigraphs. Recall that a convenient 
way (at least for theoretical purposes) to generate a random regular graph 
is the so-called configuration model or pairing model, see e.g. [1] or [19]: We 
start with nd points partitioned into n cells of d points each. We then take 
a random pairing of the points into nd/2 pairs (assuming nd to be even). 
Collapsing each cell to a vertex and regarding each pair as an edge, we obtain 
a random <i-regular multigraph that may contain loops and multiple edges; 
we denote this random multigraph by G*(n, d). (The points themselves are 
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called half-edges.) It is well-known, and easily seen, that if we condition 
G*(n,d) on being a simple graph (no loops nor multiple edges), then we 
obtain the uniformly distributed random regular graph G(n,d). Moreover, 
it is well-known that for fixed d, the probability P(G*(n, d) is simple) tends 
to a non-zero limit as n — ► oo; hence, every property that G*(n, d) has whp, 
is whp enjoyed by G(n, d) too. In particular, we may (and will) prove 
Theorem 1.1 by proving the following extension of it; we define G*(n, 2q) by 
analogy with G c (n, 2q), by choosing uniformly at random a colouring of the 
edges with n colours with q edges of each colour. 

Theorem 2.1. Theorem 1.1 holds for the randomly coloured random regular 
multigraph G* c {n,2q) too. 

We find it useful to introduce an associated bipartite graph. (This is 
really a multigraph too, since it may have multiple edges.) Given the ran- 
domly coloured multigraph G*(n,2q), add a new vertex on each edge. We 
give each new vertex the colour of the edge it bisects, leaving the original 
vertices uncoloured. Finally, we combine the q vertices of each colour into 
a single coloured vertex of that colour. This gives us a 2q-regular bipartite 
(multi)graph with n + n vertices; the original n vertices form one side of the 
bipartition, and the n coloured vertices the other. Moreover, each coloured 
vertex comes with a pairing of the edges (or half-edges) attached to it; this 
pairing shows which pairs of edges correspond to edges in the multigraph. 
We may think of the coloured vertices as having 2q attached half-edges ar- 
ranged in a circle, with each half-edge matched to the opposite one. There 
is then a one-to-one correspondence between walks in the multigraph and 
walks in the bipartite graph (of twice the length, and beginning and ending 
at blank vertices) that pass 'straight ahead' between matched half-edges at 
each coloured vertex. In particular, rainbow Hamilton cycles in the multi- 
graph correspond to Hamilton cycles in the bipartite graph that obey this 
traffic rule. 

The colours are no longer important in the bipartite graph, but it will be 
convenient to refer to the two sets of vertices in the bipartition as 'plain' 
and 'coloured'. 

Conversely, we may start with the bipartite graph, with given traffic rules, 
and obtain the original multigraph by combining the edges two by two at 
the coloured vertices. Note that choosing the bipartite (multi) graph at ran- 
dom using the configuration model (in its bipartite version, and with traffic 
rules as above given in each coloured cell) gives back the random coloured 
multigraph G*(n,2q) with the right distribution. 

We let B*(n,n;2q) denote this random bipartite multigraph with traffic 
rules as above, that is, with n + n vertices of degree 2q and with a random 
pairing of the half-edges at each vertex in the coloured part of the biparti- 
tion. (Although viewing it as a multigraph when referring to cycles etc., all 
computations are done with the equivalent configuration model.) To prove 
Theorems 1.1 and 2.1, it is enough to prove the following. 
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Theorem 2.2. The random bipartite multigraph B*(n,n;2q) whp has a 
rainbow Hamilton cycle obeying the traffic rules if q > 4, and not if q < 3. 

Remark 2.3. One may study random (regular) graphs with other traffic 
rules at the vertices. In general, we may equip each vertex of degree d with 
a (possibly directed, and possibly random) connection graph with d vertices 
representing the incident edges; the edges in the connection graph show the 
allowed connections between incoming and outgoing edges. In our case, the 
connection graph is the complete graph (no restrictions) for one side of the 
bipartition, and a matching with d/2 edges for the other side. We do not 
know of any general study, but a few examples of this type have appeared 
in the literature: 

Garmo [8, 9] studied random railways; these are regular (typically cubic) 
graphs where the vertices (representing switches) have connection graphs 
that are stars. In [10], this was extended to graphs where a random subset 
of the vertices have a star as connection graph and the rest the complete 
graph. 

Gamburd [7] studied long cycles in random oriented cubic graphs; here 
the connection graph is a directed 3-cycle at each vertex. 

3. Small subgraphs 

The small subgraph conditioning method introduced by Robinson and 
Wormald [16, 17] has been successfully applied to several problems, in par- 
ticular in the theory of random regular graphs, see e.g. [13, Chapter 9], [19] 
and [11]. (For applications to random hypergraphs, see [5, 4].) 

As often pointed out by Alan Frieze, see [4, 5, 6], the method can be 
regarded as an analysis of variance. The main idea is that we consider some 
random variable, Y say, that counts occurrences of some structure, and let 
a parameter n — ► oo. Typically, it is easy to prove that the expectation KY 
tends to infinity, but we want to show that ¥(Y > 0) — > 1. If the vari- 
ance Var(Y) is o(EY~) 2 , then the second moment method (i.e. Chebyshev's 
inequality) immediately shows the desired result. The small subgraph con- 
ditioning method applies to cases where the variance Var(V) is of the same 
order as (E Y~) 2 , by showing that the variance can be explained, up to a fac- 
tor 1 — o(l), by the interaction between the numbers of some small subgraphs 
and the random variable Y. The desired conclusion Y > whp then follows 
by conditioning on the numbers of these small subgraphs and using Cheby- 
shev's inequality on the conditioned variables. For details, see [13, Theorem 
9.12-Remark 9.18] and [19, Theorem 4.1]. We state the results there in the 
following form (an immediate consequence of [19, Corollary 4.2]). We use 
[x] m := x(x — 1) • • • (x — m + 1) to denote falling factorials. 

Theorem 3.1. Let Aj > and 5i > —1 be real numbers for i = 1, 2, . . . and 
suppose that for each n there are random variables = Xj(n), i = 1, 2, . . . 
and Y = Y{n), all defined on the same probability space Q = Q n such that 
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Xi is nonnegative integer valued, Y is nonnegative and EY > (for n 
sufficiently large). Suppose furthermore that 

(i) For each k > 1, the variables X\, . . . ,X k are asymptotically inde- 
pendent Poisson random variables with EAQ — ► Xi, 

(ii) if m = Aj(l + Si), then 



for every finite sequence mi, . . . , of nonnegative integers, 

(iii) Hi x i$i 2 < 00 > 

(iv) Ey 2 /(Ey) 2 <exp(^.A i (5 i 2 ) + o(l) asn^co. 

Then, if £ is the event A$ i= -i{Xi = 0}, F(Y > | £ ) — ► 1. In particular, if 
5 > — 1 for every i, then Y > whp. 

We will actually use Theorem 3.1 with a doubly indexed sequence X^-; 
obviously, this is just a matter of notation. 

In many applications of the small subgraph conditioning method, the 
variables Xi are the numbers of cycles of different lengths. This has perhaps 
misled some into the belief that the short cycles are expected to provide 
the answer in all cases. But they play the central role for most problems 
only because they are the only possible 'unusual' small subgraphs. The 
subgraphs of fixed size in a random <i-regular graph are very well behaved. 
Near a random vertex, such a graph looks locally like a tree. But even that 
statement can be misleading when we consider what comes shortly. The 
thing to focus on is that, because whp no two short cycles are near each 
other, the number of subgraphs of any particular type are determined by 
the numbers of short cycles. 

In our case, it will turn out that conditioning on the numbers of small 
cycles in G*(n, 2q) does not explain all of the variance of the number of rain- 
bow Hamilton cycles; we have to consider also colourings. Note that for each 
fixed length i, there are only a few cycles of length i (the expected number 
is O(l)), and whp they are all rainbow, so we would expect no explanation 
of variance to be caused by the numbers of intrinsically differently coloured 
short cycles. However, we may consider, for example, the number of paths 
of length i where the first and last edges have the same colour. The expected 
number is 0(1). These are intrinsically different from short rainbow paths, 
and it turns out that these structures too will be significant in the analysis 
of variance. 

Perhaps surprisingly, there is even more to consider. The existence of two 
short paths, each joining a blue edge to a red edge, is significant, even though 
they are two different blue edges and two different red edges in distant parts 
of the graph. However, this is not surprising given the discussion above 
about small subgraphs. Colours are clearly relevant in our present problem, 
so we should consider coloured subgraphs. Typical small subgraphs, not 
necessarily neighbourhoods of vertices, are forests with distinctly coloured 



E(Y[X 1 ] mi ---[X k ] mk ) 




(3.1) 



i=l 
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edges. The numbers of small forests in which some of the edges are coloured 
the same thus qualify as special small subgraphs. It was for this reason 
that the method was called small subgraph conditioning in [19], rather than 
short cycle conditioning. This is indeed the first application of the method 
in which the small subgraphs involved are disconnected. 

To describe the general situation precisely, we work with the random 
bipartite multigraph B*(n,n;2q) defined in Section 2, and let Y be the 
number of Hamilton cycles in the multigraph that obey the traffic rules. 
Recall that Y equals the number of rainbow Hamilton cycles in G*(n, 2q). 
Further, for each i > 1 and j with < j < i, we let X^ be the number 
of cycles of length 2i in B*(n,n;2q) that violate the traffic rules at exactly 
coloured j vertices. (Thus, Y = X n o, but we are mainly interested in for 
small i.) Note that Xio equals the number of rainbow i-cycles in G*(n, 2q), 
and thus it whp equals the number of i-cycles in G*(n, 2q), while Xn whp 
equals the number of paths of length i + 1 where the first and last edges have 
the same colour. (This holds only whp, since for Xn the endpoints of the 
path may coincide with each other or with some interior point.) We may 
similarly interpret Xy for j > 2, at least whp, as the number of certain 
collections of j paths, generalising the example mentioned above, but we 
leave the details to the reader. 

We state three lemmas that will be proven in the following sections. 

Lemma 3.2. Suppose that d = 2q > 4. Then 



Hence, asn^oo, E(Y) -> for d < 6 but E(Y) -> oo for d > 8. 

Lemma 3.3. Conditions (i) and (ii) in Theorem 3.1 are satisfied for the 
variables (Xij)ij and 



Proof of Theorem 2.2. First note that if d < 6, then E(Y) — > by Lemma 3.2, 
and thus P(Y > 0) — ► 0, i.e. Y = whp. In other words, there is then whp 
no rainbow Hamilton cycle in G*(n,2q). 

For the remainder of the proof, assume that d = 2q > 8. By Lemma 3.2, 
E(Y) — ► oo. We want to show that Y > whp. We employ Theorem 3.1 
with X^ as defined above, and Ay and Sij as given in Lemma 3.3. Note that 





Lemma 3.4. Suppose that d > 4. Then 
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Sij > — 1 for all i and j, so it remains only to show that the assumptions 
(i)-(iv) in Theorem 3.1 hold. For (i) and (ii), this is Lemma 3.3. 
For (hi) and (iv) we split the sum into two parts. 

J2 Xi i 5 l =Y.2i (j) (d-l)*(d-2)i 
i,j>0 i,j>l v/ v J y > 

1=1 

1^1// d + 2 y _L \ 

- 2 ^ i\\(d-l)(d-2)J (d-1)*/ 



1 / d + 2 \ 1, 
-In 1 - -- + - In 1 



and 



2 V (d- l)(d- 2)J 2 V d-1 

1 d 2 -4d 1 d-2 
= In — — H — In 

2 - - 2) 2 d-1 

= lln^! 
2 d(d - 4) 

VA,o4= T -7T^ = -lnfl-T^-Hlnfl + T^- 
^ 10 ^2i(d-lV V d-1/ V d-1 

j j odd 

, d-2 , d d 
= — In — h In — = In ■ 



d-1 d-1 d-2 
Consequently, 

(d-2) 2 , d 1, d 

In 



v- x x2 1 , ( d - 2 ) , , 



2 d(d-4) d-2 2 d-4 

hi 

This proves (hi), and together with Lemma 3.4 also (iv). □ 

4. Expectation 

Proof of Lemma 3.2. There are n\ 2 /2n ways to arrange the 2n vertices in 
a cycle, with plain and coloured vertices alternating, and for each such 
arrangement d(d— 1) ways to choose the half-edges at each plain vertex and 
d ways to choose the half-edges at each coloured vertex (obeying the traffic 
rules). For each such choice, the probability that the selected An half-edges 
are connected to each other in the specified order equals ((d — 2)n)\/(dn)\. 
Consequently, using Stirling's formula, 

d 2 "(d-l)"n! 2 ((d-2)n)! _ ^ d 2 (d - l)(d - 2)^ 2 y _ 

2n(dn)\ ~ U \ d d ) - U y JW > ' 

where /(d) := (d- 1)(1 - 2/d) d - 2 . We have /(4) = 3/4 < 1, /(6) = 80/81 < 
1, /(8) = 5103/4096 > 1, and /(d) > (d - l)e~ 2 > 1 for d > 8. □ 



8 



SVANTE JANSON AND NICHOLAS WORMALD 



5. Short cycles 

Proof of Lemma 3.3. We use arguments that have become standard for sim- 
ilar problems for random ragular graphs, see e.g. [13, Section 9.4] or [19, 
Section 4.2]; we will thus omit some details. 

For (i), we use the method of moments. It suffices to show that 

ij i=l 

for every finite set of non-negative integers {to^}. For convenience, we will 
only treat the expectation of a single as in all similar problems, the 
argument extends immediately to (mixed) higher factorial moments. 

To calculate EJQj, we count the appropriate oriented cycles with a des- 
ignated initial vertex, which we require to be plain; this counts each cycle 
2i times. The vertices in the cycle may now be chosen in [n]| ~ n 21 ways. 
Consider first the case j = 0, i.e. cycles obeying the traffic rules everywhere. 
For each choice of vertices there are, as in Section 4, d{d — 1) ways to choose 
the half-edges at each of the i plain vertices and d ways to choose the half- 
edges at each of the i coloured vertex. Finally, the probability of pairing the 
4i chosen half-edges into 2i edges is l/[dn]2i ~ (dn)~ 2t . Hence, 

2 iE ^ = d2, <<;- 1 i > , w. 2 ^ p.^. 

[dn\ 2 i 

For j > we argue similarly. The traffic rules are to be violated at precisely 
j coloured vertices. These may be chosen in Q) ways, and at each of them 
there is additional factor of d — 2 for the choice of the out-going half-edge. 
Hence we obtain, for all i and j, 

or E Xij — ► Ajj . 

For (ii), we first observe that the left hand side of (3.1), by symmetry, 
remains the same if we fix two half-edges at each vertex, always choosing 
two opposite half-edges at the coloured vertices, and then replace Y by the 
indicator that the chosen half-edges comprise a rainbow Hamilton cycle. 
Denoting this event by Hi, we thus want to show 

k 

-(II -v,, - |«i) (s-i) 

ij i=l 

Condition on Hi, and let H\ be the (unique) rainbow Hamilton cycle 
that uses the chosen half-edges. It is easily seen that the remainder of the 
graph G*(n,2q) can be regarded as the random multigraph G*(n,2q — 2), 
and that this is independent of H\. Hence, the left hand side of (5.1) equals 
the expectation in the union of a random rainbow Hamilton cycle H\ and 
an independent G*(n, 2q — 2) on the same vertex set. 
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For the same reasons as in (i), we will only consider a single expectation 
E(Xij | Hi). Consider first the case j = 0. We may, as for (i), choose the 
vertices of the 2i-cycle in n 2i (l + o(l)) ways. We then decide whether the 2i 
edges are in the Hamilton cycle Hi or in G*(n, 2q — 2); we denote the choices 
by a s G {1, 2} for s = 1, . . . , 2i. At a plain vertex where the incoming edge 
is to have type a G {1,2} and the outgoing edge type /3, there is for each 
possible incoming half-edge a Q/ g choices of the outgoing, where the numbers 
a a p are conveniently collected in the matrix 



A:=(a a p)=(l d d _fj 



At the coloured vertices there is only one choice for the outgoing edge, and 
it has to have the same type as the incoming; we encode this as b a s with 
B := {b a p) = (q i), the identity matrix. Note that we have not counted the 
number of incoming half-edges; this is because these numbers cancel when we 
take into account the probability of making the connections; the probability 
of connecting a half-edge of either type to some incoming half-edge of the 
same type at a given vertex (in the opposite part of the bipartition) is 
(l + o(l))n™ 1 . Moreover, the probability that all 2i connections are made is 
(1 + o(l))n -2 *, except in the case when all a s = 1, which is impossible for 
n > i. Consequently, 

2 

" Oll&ll 

ai,-,«2s=l 

= Tr(ABy - 1 = Tr(^) - 1 = (d - 1)* + (-1)* - 1, 

since A has the eigenvalues d—1 and —1. We have thus shown (5.1) for this 
case, with 

_ (d-i)t + (-i)i-i 

ViO — ^ — AjoU + Oio) 

as required. 

For j > we argue similarly. Now we have to choose j coloured vertices 
where the traffic rule is violated, and for these vertices the matrix B is 
replaced by 

'0 d-2^ 



B:=A-B 



2 d-4 



Luckily, A and B = A—I commute, so all Q) choices of the violating vertices 
give the same result, and thus, for j > 0, 

2iE(X tj | Hi) Q TriAtW) = Q Tr(A*(A - 

= Q((d-iy(d-2y + (-iy(-2y), 

which equals 2iA^ (1 + 5ij) as required in this case too. □ 
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6. Variance 

In this section we prove Lemma 3.4, thus completing the proof of Theo- 
rem 1.1. 

We compute E Y 2 by calculating the probability that a given ordered pair 
of Hamilton cycles (H 1,-^2) are contained in the pairing corresponding to 
B*(n,n;2q), and summing over all possible ordered (H\,H2). 

The first part is similar to the treatment of Hi in Section 5. Let k denote 
the number of coloured vertices in which the same half-edges are used by 
both Hi and H2, and let j denote the number of blank vertices of this type. 
If Hi / H2 then the half-edges shared by the two cycles occur in k — j 
"strings" of consecutive half-edges around H\. (Note that each string ends 
at plain vertices.) The strings also occur in H2, though in a different order. 
Together, Hi and H2 determine the pairs containing two of the half-edges at 
k coloured vertices and four at all other coloured vertices, so 4n — 2k pairs 
in all (as each pair contains just one coloured vertex). Hence, defining 7ii 
as the event that Hi occurs, 

P(W2 | Wl) .««-4>" + ^ 



((d-2)„)\ ' 
and we may write 

EY 2 _l + Y, k , J<n N(k,j)nH2\Hi) 



(EY) 2 EY 

where N(k,j) is the number of different H2 (or, to be precise, sets of pairs 
corrsponding to H2) overlapping any given cycle Hi with particular values 
of k and j. The term 1 accounts for the case H2 = Hi, when k = j = n. 

Since EY was given precisely in Section 4, all that is left to evaluate is 
N(k,j). Note that the cardinality of the set Si of plain vertices with three 
half-edges contained in H1UH2 is 2k — 2j, and for the set S2 of plain vertices 
with four half-edges, it is n — 2k + j. 

By elementary counting, the number of ways to place k—j strings referred 
to above onto Hi is 

n fk\ fn — k — 1 N 
k \j)\k-j-l, 

Here the first factor converts the problem from a cyclical one to a linear one 
in which the first coloured vertex is one of the k special ones. The second 
factor is for deciding the relative positions of the k coloured vertices in a 
sequence of k — j strings and the third is for deciding the positions of the 
vertices not in strings. 

The number of ways to choose the half-edges being used by H2 at each 
plain vertex is (d - 2)^((d - 2)(d - 3))l 52 l, and (d - 2) n ~ fc for the half- 
edges at the set S3 of coloured vertices not used by Hi (obeying the traffic 
rules). These choices determine the direction that H2 passes through each 
vertex except for those in Si , which are determined by the direction it passes 
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through each string. These directions can be chosen in 2 k J ways, giving 



ways to make these choices. (For convenience we will count oriented versions 
of H2 and will divide by 2 at the end.) 

Now that the order of half-edges used by H2 is determined at each vertex, 
it remains to choose the remaining 2(n — k) pairs. These pairs must connect 
the strings and the vertices in S2 U S3 into a Hamilton cycle in a bipartite 
fashion (regarding each string as a vertex) and connecting the 'out' half-edge 
at a vertex to the 'in' one at the next. The number of such choices of pairs 
is 

(n - k)\(n - k - 

Multiplying all the displayed factors together and dividing by 2 to un-orient 
H2 gives N(k,j). Combining with the earlier equations then produces 



f(n,d,k,j) = 

n 2 (k - 1)! (n - k)l 3 2 k -i(d - 2) 2n ~ k ~i{d - 3) n ~ 2k+ ^ ((d - 4)n + 2k)\{dn)\ 
(n - k) 2 (k - jj\ (k-j- 1)! jl (n-2k + j)\ (Jd - 2)n)\ 2 d 2n {d - l) n nl 2 

As usual, upon applying Stirling's formula we find that the powers of n/e 
cancel and we are left with 



2 k - j (d - 2) 2n ~ k - j (d - 3) n ~ 2k+j 




k<n j<k 



(6.1) 



where 




(6.2) 



where 



fo(n,d,k,j) = 



^(dn - 4n + 2k)nd 



2n(d - 2)^k{n - k)j(n - 2k + j) 




(6.3) 



T = {((5, a) : < a < 5, 5 + a < l}. 



The partial derivatives of In F are 

d(lnF) t 2 (t-l)5 3 (l-5-a) 
85 ~ (t- 25) 2 (1 -5)(S- a)' 
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g^ = h »('-«)(i-«-«), (M) 

aa (i — l)ar 

Setting these equal to zero gives necessary conditions for a stationary point 
of F: 

t 2 {t - 1)5 3 {1 -5-a) = (t- 25) 2 (1 - 5)(5 - a) (6.5) 

and 

(t + l)a 2 - 2ta + 2i<5(l - 5) = 0. (6.6) 
Solving the first equation for a and substituting this into the second shows 
that the value of 5 at a stationary point in the interior of T must be a root 
of 

g (S) = (1 - 5)((t - 25) 2 - t 2 (t - l)5 2 f - 2t 3 5(l - 25) 2 (t - 25) 2 , (6.7) 
which is quintic in 5. Now g factorises as 

g (5) = (t + 2)(5 -5)h(5) (6.8) 

where 

h(5) = t 5 5 4 - 25{25 3 + 5 2 - 25 + l)t 4 + (95 4 - 125 3 + 65 2 + l)t 3 

+ 2<5(<5 - 1){35 2 + 25 + 3)t 2 - 45 2 (5 + 3)(<5 - l)t + 85 3 {5 - 1), 

and So = t/(t + 2) which, as we will show, determines the maximum value 
of F in T. 

The second derivative of h with respect to 5 is 

12t 5 5 2 + (-A85 2 - 125 + 8)t 4 + 12(9<5 2 -65 + l)t 3 

+ (725 2 - 125 + A)t 2 + (-A85 2 - 485 + 24)t + 965 2 - 485. 

This can be rearranged as 

I2(t - 6)t 4 5 2 + (245 2 - 125 + 8)t 4 + 12(9<5 2 -65 + l)t 3 + 125{t 2 - 4) 

+ 125 2 (t 2 - 6t) + 4(15<5 2 -65 + l)t 2 + 24(5 2 - 25 + l)t + V65 2 , 

in which each collected term is clearly nonnegative for t > 6 and all 5 > 0. 
Thus h is a convex function of 5 for each t. 

From (6.7) we compute firstly 5(1/2) = (t - l) 2 (t - 2) 4 /32 > 0, secondly 

9(1/ Vi) <t 2 - 2f 5 / 2 (l - 2/Vi) 2 {t - l) 2 < 

since t > 6 and 25(1 - 2/ v / 6) 2 > 0.8, and thirdly, since t > 6 implies 
t 2 - 7t + 8 > and thus l-l/t< 2(1 - 2/t) 2 , 

g(l/t) < (1 - l/t){(t - 2/t) 4 - t 2 (t - 2/t) 2 ) < 0. 

Note that h has the same sign as g for < 5 < 1/2 as <5o > 1/2. Thus h(l/2) 
is positive and h(l/y/t) and h(l/t) are negative. So by the convexity of h, 
the only zeros of h for < 5 < 1 lie in < 5 < 1/t or 1/y/t < 5 < 1/2. We 
show separately that these two subsets of T can hold no stationary points 
of F. 

Case 1: < 5 < 1/t 
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Substituting a = (35 into In F, and taking the second derivative with respect 
to 5, we obtain 

d 2 In F((35, 6) _ 2t - (3t - 4t5 + 2t5 2 - 2(3td + 2(3t5 2 + 2(35 
d5 2 ~ (1 - 5)(1 - 5 - (35)5(t - 25) 

Since a < 5 in T, we have (3 < 1 and so the factors in the denominator are 
all positive. The numerator is at least 2t — t — 4t5 — 2t5 >0as5<l/t<l/6. 
Hence, In F can have no local maximum in T for such 5. 

Case 2: < 5 < 1/2 

For such 5, from (6.6) we obtain 2ta > 2t5(l — 5), i.e. a > 5(1 — 5). So, 
using (6.5), at a stationary point 

(1-5) 2 _ 1-5-5(1-5) 1-5 -a _ (t - 25) 2 (1 - 5) 
I 2 ~ 5-5(1-5) K 5-a ~ t 2 (t - 1)<5 3 

and so, since 5 2 > 1/t, 

1-5 (t - 25) 2 
~1~ < t(t- 1) ' 

But this fails at 5 = 1/2, and the derivative of the left hand side with respect 
to 5 is less than —4 for 5 < 1/2, whilst that of the right is easily greater 
than —4. So the inequality fails, and there is no such stationary point. 

We conclude that 5 = 5o determines the unique local maximum in the 
interior of T. The boundary of T must also be investigated. Consider- 
ing (6.4), there is no local maximum at a boundary point with < 5 < 1, 
since d(ln F)/da tends to oo as a tends to from the right, and to — oo 
as a tends to min(5, 1 — 5) from the left. A similar argument applies to 
eliminate (a, 5) = (0, 1) from consideration, since moving along the bound- 
ary where a = 1 — 5, F is a smooth function times a~ a . This leaves 
only the point a = 5 = 0, which is indeed a local maximum, with value 
F(0, 0) = (1 + 2/t) t /(t + 1) < 1 for t > 6. 

To deduce that 5o determines the unique global maximum in T, we only 
need to observe that the corresponding value of a is «o = 25o/(t + 1), and 
that F(ao, 5o) = 1. 

The rest of the argument is totally standard for such variance calcula- 
tions, as in [6] for example, so we omit the justifications. The point (ao,5o) 
corresponds to k = kq = 2/d, 7 = Ko/(d — 1). Putting k = kq + k/^/n and 
7 = To + t/v^' an d expanding ln(G") (G defined in (6.3)) about k = 7 = 0, 
we find up to quadratic terms in k and 7 

G n ~ e Cl ^ 2+C2 ^ +C3 ^ 2 
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where 

d(d 3 - 3d 2 + Ad + 4) 



ci 



C2 



C3 



4(d-2) 2 (d-3) 
d(d- l) 2 



(d-2)(d-3) 
d(d- l) 2 



4(d-3) ' 

Here C3 is clearly negative, and the determinant Z) = 4C1C3— c 2 of the Hessian 
of the quadratic form is positive, as we expect since the expansion is at a 
local maximum. The routine argument now gives from (6.2) 

^2 /( n ' rf ' fc 'i) ~ fo{n,d,nKo,nj )2im/y/D 

k,j<n 



d 



d-4 

Recalling (6.1) and that EY — > 00 now establishes Lemma 3.4. 

7. Rainbow matchings 

In this section we briefly consider the analoguous problem of the existence 
of a rainbow perfect matching in a randomly coloured random regular graph. 
We will omit the details of the calculations. 

The model is now slightly different. We consider a random regular graph 
G(2n,d) with an even number of vertices (d may now be arbitrary), and 
colour randomly the nd edges with n colours, d edges of each colour. We 
then ask whether there exists a rainbow matching consisting of n disjoint 
edges of different colours. 

We can translate this to a random bipartite (multi) graph as above; now 
the bipartite graph has 2n plain vertices of degree d and n coloured vertices 
of degree 2d. Let Z be the number of rainbow perfect matchings; in the 
bipartite version, Z is the number of decompositions of the graph into n 
disjoint paths of length 2, with 2 plain and 1 coloured vertex each. 

Calculations as above yield 

EZ= (2n^(2rf-2^ = ( n */V-Vr->y\ 
(2dn)! V ^ d ' J 

and it is easily checked that, as n — > 00, E Z — > for d < 6, while E Z — ► 00 
for d > 7. In particular, for d < 6 there is whp no rainbow perfect matching. 

Furthermore, for d > 7, an argument similar to the one in [2] (and much 
simpler than the proof of Lemma 3.4 above, since we only need to maximize 
over one variable) yields 

E(Z 2 ) _ d- 1 
WW ~* 7d(d^3) ' 
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Finally, denning as before, Theorem 3.1 applies when d > 7 with 



Hence there exists a rainbow perfect matching whp when d>7. 

By analogy with the open problem in Remark 1.2 one might further 
ask whether there exists a decomposition into d rainbow perfect match- 
ings. Computing the expected number of such decompositions reveals that 
when d < 11, whp no such decomposition exists. This compares with the 
corresponding result for uncoloured graphs, that G(2n, d) whp has a de- 
composition into d perfect matchings as soon as d > 3 [12, 15]. 
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